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Raman cooling of non-zero-spin atoms in the presence of gravitational and external magnetic 
fields is investigated. The magnetic field is adjusted so as to compensate for the gravitational force 
acting on ground-state atoms. The dark state (DS) is created and supported in momentum space 
with additional velocity-selective two-photon transitions. The minimum allowed temperature is 
found to be determined only by the width of velocity selection and therefore can be much less than 
the gravitational limit. A complete set of analytical formulas describing cooling of a dilute atomic 
sample is derived. They serve as the basis for numerical simulations which are carried out in the 
one-dimensional (ID) case. 
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I. INTRODUCTION 

Methods of laser cooling of atoms have progressed dramatically in recent years. There are three typical temperature 
scales characterizing various methods. The first one is specific to the most common scheme in which Dopplcr shift 
causes the radiation pressure force to be velocity dependent, thus damping atomic motion when the laser frequency 
is tuned below an atomic resonance. The minimum temperature Try for atoms cooled in such a way is known as the 
"Dopplcr limit" . It is proportional to the natural width of the laser-driven transition 7 ksTo = ^7/2, and for the 
D line of Na the Dopplcr limit is approximately 240 /iK. 

Schemes based on dissipation of atomic energy via interaction with the vacuum modes of electromagnetic field have 
a lower limit on the achievable temperature defined by the minimum of recoil energy which an atom obtains after 
spontaneous photon emission. The corresponding scale is known as the recoil limit £;bTr = (hk) 2 /2M , where k is the 
wave number of emitted light, and M is the atomic mass. For Na it approximately equals to 1 fiK. 

To overcome this limit two subrecoil cooling methods have been developed and demonstrated: velocity selective 
coherent population trapping (VSCPT) || and Raman cooling || . Both methods imply the existence of the so-called 
dark state, which does not interact with light, has a long lifetime and occupies only a few modes in momentum space 
H . During the cooling cycle atoms diffuse into this state due to random recoil induced by spontaneous emission and 
accumulate in it. Since DS has a vanishing absorption rate of light, the final temperature is restricted not by the 
recoil limit but by the time of cooling which, however, cannot be greater than the lifetime of DS. 

In the absence of Earth gravity, infinitely long cooling times would be possible. In practice the gravitational 
field pushes atoms from DS, reducing its lifetime dramatically. For any 3D configuration this defines the third 
characteristic temperature scale || ksTc = Mg/k, the gravitational limit. It lies below the recoil limit for most 
atoms, e.g., T G = 0.07Tr for Rb, and T G = 0.003T fl for Na. 

Two ways are envisioned to prepare a stable quantum state of matter in the gravitational field: to bound particles 
or to suspend them free in an inhomogeneous magnetic field using Stern-Gerlach effect. In the first approach atoms 
are confined by a conservative trapping potential which can be realized, e.g., in a far-off-resonance or a dipole trap, 
where an intensity gradient provides a spatially dependent ac Stark shift. In momentum space, up to now only the 
existence of an approximate dark state has been demonstrated ||, characterized by a decay rate in a special ID atomic 
and laser field configuration much smaller than that of all other states in the trap. The finite lifetime of approximate 
DS evidently restricts the cooling possibilities in a trap, leaving the question about going below the gravitational limit 
to be clarified. However, a scheme which is based on the creation of a dark state in position space with the help of 
an appropriate spatial profile of the cooling laser, e.g., in a doughnut mode, seems to be much more efficient, allowing 
to cool a significant fraction of atoms to the ground state of the trapping potential. 

Another approach may be applied to atoms possessing a magnetic moment. Superimposing a weakly inhomogeneous 
magnetic field onto the path of pre-polarized particles and appropriately adjusting the field gradient it is possible to 
compensate the effects of gravity for a definite internal atomic state. However, the magnetic field induces spatially 
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dependent shifts of the Zeeman levels, which lead to unwanted residual excitation from the DS in the framework of 
any traditional subrecoil cooling method. Moreover, in the case of VSCPT the dark state cannot be an eigenvector of 
the total Hamiltonian since only one of the internal states forming the superposition which is not coupled to the laser 
field may escape gravity. Thus, both VSCPT and Raman cooling mechanisms in their standard form are incompatible 
with the last approach. 

To resolve this problem we suggest a modification of Raman cooling method, in which the ground-level atoms are 
made motionally free with the Stern-Gerlach effect and the DS is created and supported in momentum space of these 
atoms with additional velocity-selective two-photon transitions. The transitions couple external momentum states of 
the same ground internal level and are organized in such a manner that DS cyclically occupies different thin sets of 
velocity modes while remaining unreachable for the Raman excitation-repumping pulse sequences at all times. 

In Sec. [n] a detailed qualitative treatment of the suggested scheme is given. For reasonable experimental conditions 
all the stages of the scheme admit analytical descriptions which are presented in Sec. Ill - Sec. |v|. Specifically, 
the formulas describing a coherent two-photon transition when an atom is placed in a superposition of two plain 
electromagnetic waves with arbitrary directions of the wave vectors are presented in Sec. [II. In Sec. |Ev] Raman 
excitation to the closest hyperfine level is investigated in the regime in which the photon spontaneous emission may 
be neglected. Both the exact quadrature and convenient approximate expressions are derived. In Sec. [v| the optical 
pumping of atoms back to the ground state is considered for short times of light pulse, i.e., when the external potential 
field does not affect the ato mic m otion substantially. A numerical simulation of 100 cooling sequences in one dimension 



is given in Sec. VI. Section VII concludes with a summary of the obtained results. 



II. QUALITATIVE TREATMENT 

Consider for definiteness an atom with a, J — ^ to J — ^ transition, e.g., sodium or cesium. The magnetic field B(r) 
applied to compensate the gravity is supposed to contain a homogeneous component Bo directed along the gravity 
acceleration Bo tt g- The remaining inhomogeneous part of the field Bi(r) ~ B(r) — Bo should be small compared 
to this component, 

|Bi(r)|<flo = |Bo|. (2-1) 

As we will see below, to fulfil this condition it is necessary to take Bo in the range 10 3 -j- 10 4 G. In practice such a 
field is strong enough to induce Zeeman shifts which considerably exceed the hyperfine splitting intervals ~ Hujufs 
(but not the multiplet ones). Therefore an internal atomic eigenstate \J,I,Mj,nij) may be well described using the 
set of quantum numbers consisting of the angular momenta of the electronic shell J and the nucleus /, and their local 
projections Mj, mj on the direction of the magnetic field. 

In the framework of perturbation theory, | J, I, Mj, mi) represents a combination of eigenstates | J, I, Mj, m/)(°) 
related to the atomic Hamiltonian without the hyperfine interaction, 

| J, I, Mj, mj) = \J,I, Mj, mi) ^ + - - 

2/j,b9lBo 

x {[(J + Mj)(J - Mj + + mj + 1)(J- mj)] 1 ' 2 

X\J,I,MJ- 1,771/ + 1) (0) 

— [(J + Mi + 1)(J- M J )] 1 / 2 [(7 + mj)(I - mi + I)] 1 ' 2 

x |J,/ > M J + l,n»j-l)W}, (2.2) 

where a is the hyperfine coupling constant (a cx Iiljhfs, e.g., for Na a/% = 885.8 MHz) and t/z, denotes the Lande 
factor. The corresponding energy eigenvalue is determined not only by the multiplet level Ej but also by the magnetic 
field B(r) = |B(r)| and therefore is spatially dependent 

E\J,i,Mj,mi)( r ) = E J + aMjmj 

+ {hb9lMj - /j nuc mi)S(r), (2.3) 

where /z nuc is the nuclear magnetic moment. Because of the condition ( |2.l| ) such a spatial dependence, however, 
mainly arises from the longitudinal (b\{y) — Bo • Bi(r)/£?o), rather than the transverse (B^r)) component of the 
vector Bi(r), provided that the components are defined relative to Bq. This is evident from the expression 
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B(r) = ^[B Q + Bl(T)} 2 +[B^(r)] 2 

~B + i?!(r)+ [B^(r)] 2 /(2So), (2.4) 

where the term containing Bj'-(r) is small and can be neglected. Consequently, by adjusting the gradient of the field 
i?['(r) one can achieve translational invariance of the ground state |1) = |l/2, I, —1/2, I) in three dimensions: 

E w (r) - Mg • r = const. (2.5) 

For example, to balance the gravitational force in this way for sodium it is necessary to create a gradient Vs| (r) = 
&ig/|g|, where b\ = —4.033 G/cm. This condition does not contradict the Maxwell equation V • Bi(r) = , be cause 



variation of B^-(r) is not restricted. Note also that the choice Bq = 10 3 -j- 10 4 G maintains the condition (2.1) very 
well within a spatial region of the size ~ 10 cm. 

All the other levels are affected by the residual external potential. In particular, after a transition from |1) to the 
neighboring state |2) = 1 1/2, I, —1/2, 7—1) the atom experiences a force 



M/j, nuc g 



(l/2)(i B g L + 



(2.6) 



In our scheme, we use pulses of laser light at frequencies oj\ and LO2 which are roughly tuned to the |1) — > |3) and 
|2) — * |3) transitions, where |3) = 1 3/2 , 1, — 3/2, 1) is an excited state with the lowest energy. The typical size 2L of 
atomic sample is restricted by the condition L <C a/(Mg), which allows to regard Ei s \(r) as the closest to resonance 
excited level within the whole i nter action domain. Indeed, the force f3 acting on the atoms in the state |3) may be 



estimated from Eqs. ( |2.3| ) and fl2.5|) as | £"3 1 ~ Mg. The maximal spatial shift of the level ~ MgL which it induces 
is much less than the hyperfine splitting intervals (MgL C a ~ huj-ups), and the hierarchy of detunings is retained. 
Therefore an atom initially in |1) or |2) state behaves as a three-level system with respect to the processes with 
stimulated emission of photons. 

Since the atomic dipole momentum operator d is diagonal in quantum numbers I and toj in the basis 
I J, I, Mj, mi)(°\ the transitions which change m/, e.g., |2) — > |3), are allowed only due to hyperfine interaction, 



as is seen from Eq. (2.2). The value of any matrix element like |(3|d|2)| is approximately cx 7jHp|(3|d|l)|, where 
r/HF = a/(2/i_B.9L-Bo) Cf As a consequence, the upper state |3) decays to the lower ones preferentially in the channel 
|3) — * |1) (with the rate 7). This circumstance makes it possible to deal with an atom as a three-level system even if 
spontaneous photon emission takes place. 

When the atom is irradiated with two laser beams at frequencies uj\ and u>%, the two-photon Raman transition from 
|1) — > 1 2) has twice the Doppler sensitivity of a single-photon transition provided that u)\ — lo- 2 ~ whfs an d the beam 



wave vectors ki, k2 are opposite ||. However, if we take into account the force (2.6), a wide set of atomic momenta 
p may satisfy the resonance condition, as follows from the energy conservation: 

TiAi - 2p • A p /M = HA 2 - f 2 • r + 2A 2 p /M. (2.7) 

Here detunings A m = uj m + [E\ m \ (0) — E\§\ (0)]/7i, m — 1,2, are defined in the center of atom-laser interaction region 
(r = 0), A p = h(ki — k 2 )/2, and A p = \A p \. The dip in the velocity dependence of absorption rate broadens so that 
the width of the trapping zone || becomes 

Sv ~ L|f 2 |/(2A P ). (2.8) 

As a consequence, since the sample of unconfined particles considered in this paper may spread up to L ~ 1 cm 
during the cooling, the effective temperature of atoms left in the state |1), which constitutes ~ M(5v) 2 /(2fcs), 
generally lies far above the gravitational limit. For example, in the case of sodium, where A p /% = 1.07 x 10 5 cm -1 
and |fa|/^ = 7.3 x 10 4 cm -1 s _1 , such a temperature may reach 6.4Tg. 

Despite insufficient velocity selectivity of the |1) — > |2) transition, state |2) may be used in Raman excitation cycle. 
To avoid unwanted radiation impact on the selected group of particles, which are referenced here as the DS atoms, 



one should move them in momentum space to another place, where the resonance condition (2.7) brakes down. It can 
be achieved by means of a two-photon |1) — > |1) transition while the atom is irradiated with two noncolinear laser 
beams at the same frequency Wi. 

If the ground- level initial momentum distribution along the direction of vector A p were as shown in Fig. 0(a) , such 
a transition would have selectively brought particles concentrated near the point — A p (the DS, as we will see below) 
to the point A p , and vice versa. To prove this imagine an atom with momentum p passing through a superposition of 
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two laser beams. The superposition may be treated as a diffraction grating in the case ki || k2 (standing wave) fl9"|.[To||, 
or as an effective atomic hologram when directions of the wave vectors are arbitrary |llf| . At low laser light intensity 
and large detuning Ai only the first-order Bragg scattering is of importance fl^ ]. In this case, two diffraction modes 
with indices and 1 resonantly couple with each other jl(J[l3| . Physically, the first-order Bragg resonance corresponds 
to an absorption and stimulated photon emission process from one laser beam to another. As a consequence of the 
atomic kinetic energy conservation one gets the Bragg resonance condition 

±p-A p = A2, (2.9) 

which is satisfied for any momentum with the component p — ±A p along the vector A p . Figure 0(b) contains the 
final distribution, the peak around A p being the moved DS. So the first step of our scheme consists in the momentum 
transfer of DS as it is indicated with arrows in Fig. 0(a). 

In the second step of cooling, the Raman excitation cycle [^) takes place. In accordance with Eq. ( |2.7| ), atoms with 
any negative p can be transferred to state |2) by varying the difference of beam frequencies. Due to the finite width 
of trapping zone atoms with positive p < MSv also have a chance to undergo transition. The DS, being hidden near 
the point p = A p > MSv, does not take part in this process, as illustrated in Fig. 0(b). 

In the third step, an optical pumping pulse at frequency 0J2 is used to return the atoms back to the state |1). It 
is important that the ground level appears to be far off resonance and laser light does not affect DS directly. The 
population of DS rises during the spontaneous emission process, which randomizes the atomic momenta [see Fig. 0(c)] . 

Then the sequence of steps 1 - 3 is repeated with opposite directions of ki and k2 involving residual positive- 
momentum atoms of the ground level in DS filling and finishing a ID cooling cycle along A p . After this stage the DS 
occupies its initial place near the point — A p . 

By choosing linearly independent vectors A™ in a set of two (m = 1, 2) or three (m = 1, 2, 3) ID cooling cycles one 
can proceed with decreasing the temperature in two or three dimensions by repeatedly applying such sets. 

To increase the efficiency of DS filling one can admit several Raman and optical pumping pulses, i.e., a number of 
steps 2 and 3, between two consecutive first steps. It can be done, for example, as in the classical method ||, where 
every Raman transition is followed by the optical repumping, or by applying a series of cycles, each including multiple 
Raman and one optical pumping pulses. 

Since the time necessary to collect all the atoms in DS is, generally speaking, infinitely long, it may be useful to 
separate the DS from background with the final first-step transitions (on one for each dimension) so that the DS 
and background atoms will move in opposite directions and eventually will not spatially overlap. In particular, when 
vectors A™, m = 1, 2, 3, form an orthogonal basis, our scheme will produce a cooled atomic beam with the average 
momentum A p + A^ + A^ as follows from Eq. ( |2.9p . The minimum allowed temperature (but not the intensity) of 
such a beam is obviously determined by the width of velocity selection specific to first-step transitions and therefore 
can be much less than the gravitational limit. 



III. GROUND-STATE TWO-PHOTON TRANSITIONS 

In contrast to the case of Bragg scattering jl0],|l2| , where diffracted modes are assumed to be spatially resolvable 
at some distance from the light standing wave, the present paper deals with short interaction times when an atom 
moves inside a superposition of two laser beams from the beginning to the end. This allows one to represent each 
beam as a plane electromagnetic wave (m = 1,2) 

E m (r, t) = E m exp(ik m • r - iu) m t) + c.c, (3.1) 

where E m stands for the complex amplitude, and 0J2 = lo\. 

To simplify the consideration the coherent scattering processes are assumed to dominate the spontaneous emission, 
i.e., the regime |Ai| ^ 7 is kept [Q. Under such a condition the one-particle density matrix jl2| has an obvious time 
evolution 

Pab(PuP2,t) = J dpi J dp 2 ^G Qa <(pi,pi,i) 

a'h' 

xG* w (p2,p' 2 ,t)p a ,v(p' 1 ,p' 2 ,t = 0), (3.2) 

where indices a,b . . . denote the internal atomic states and G aa ' (pi, Pi, t) is the Green function of the two-component 
Shrodinger equation describing atomic dynamics during the |1) <-* |3) transitions. 

In the rotating wave approximation the equation for slowly varying in time ground- and excited-level wave functions 
ipi(p,t) and ip3 (p,t) takes the form 
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[t(p) + A x ]^i(p, t) - niMp + nki, t) 
-n;v>3(p + frk 2 ,t), 



(3.3a) 



[t(p) - if 3 • V]^ 3 (P, t) - SIMP - ^1. *) 
-02*01 (P- ftk 2 ,t), 



(3.3b) 



where fi m = (3|d ■ E TO |l)/7i, m = 1,2, are the Rabi frequencies, and the terms i(p) = p 2 /(2Mh) and — zf 3 ■ V arise 
in momentum space from the kinetic and potential energy (—fa ■ r) correspondingly. 

For the situation at hand, the upper electronic state can be adiabatically eliminated from Eqs. ( |3.3a ), ( 3.3b| ) 
provided that the detuning Ai is large enough [§l|[l5 16 



|Ai| » l^xl, \tl 2 \,\f 3 \L/K. 



(3.4) 



The route by w hich one can do it implies a self-consistent assumption \ip$\ <C leading to the zero-orde r solu tion 
oftheEq. (3.3a): tpi(p,t) ~ exp{— i[t(p) + A.±]t}ipi(p, t — 0). After substitution of this expression into Eq. (3.3b) the 
latter may be solved in the framework of perturbation theory developed with respect to the potential energy term. 
In this case, the excited-level wave function acquires a representation 



t{p) - t( P - hkx) - Ai 
t(p) - t(p - hk 2 ) - Ai 



(3.5) 



where the dots denote omitted terms which include a small (oc |fs|.L/|7iAi|) first-order correction to ip3(p,t) and also 
summands which oscillate with the non-resonant freq uenc y t(p) and therefore give a negligible contribution when one 
uses the above expression within the context of Eq. ( ]3.3a| ) 



For an ultracold atomic sample one can further neglect the kinetic energy terms in the denominators of Eq. (3.5) 
so that after introducing of a new set of functions (n £ Z) 

^ n \p,t) = MP + (2n - l)A p ,t), (3.6) 

Eq. ( |3.3a ) becomes equivalent to an infinite system of equations defined in the domain T> = {p : \p ■ A p | < A^}: 

. d 



i -^ (n) (p,t) = t (n) (p)V ( " ) (P,i) 



+<^ (n+1 Hp,*)+sV (n_1) (p,f). 



where 



tW(p)=t(p+(2n-l)A p )+Ai 



Ax 



Ai 



(3.7) 



(3.8) 



and g = £l*£l2/Ai stands for the effective Rabi frequency. 

At lo w |g I only the two functions with n = and n = 1 have a possibility to influence each other resonantly in the 
system ( p/7| ) because only t^°\p) and t^\p) may be equal when p £ V. If we take into account the coupling of other 
functions, all ^ n \p,t) will get corresponding correction s oc 8^ = \g\/(t^ n \p) — t^ m \p)), where m — nil, and 



rn + n =/= I. Therefore it is possible to truncate relations (3.7), having in mind that 5^ <C 1 when the effective Rabi 
frequency is small enough. In this case, the equations for ip^(p, t) with n = 0, 1 become homomorphic with the rate 
equations describing a two-level atom, and their solution is well known (see, e.g., |l^| ). The remaining non-resonance 
functions simply undergo a free evolution. 

However, as a general rule, the original wave function ^>i(p, t) reconstructed in accordance with the formula ( ]3.6| ) 
appears to be discontinuous along the planes p • A p = 0, i2A p . To recover a smooth behavior, one can modify the 
reconstruction prescription, e.g., 



MP, t) = ^> (p + A p ,t)+ ^ (p - Ap, t) 
- exp[-it< >(p + A p )#i(p,t = 0), 



(3.9) 
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where the solutions if>^ 0,1 '{p, t) of the truncated system (3.7) must be analytically continued into the whole momentum 
space. It is easy to check that ^i(p, t) built in such a way coincides with the exact representation via ^"^(p , t) up t o 
an error of order 8^ ~ 8^ ~ \g\/t(2A p ) and, consequently, obeys (with the same accuracy) the Eqs. ( |3.3a[ ) - ( |3.3b| ). 
As a result the ground-state component of the Green function is given by 



Gii(p,p',t)= E ^ a) (P,t)S 3 (p + 2aA p -p') 

a=0,±l 

x exp[iujxt - (i/K)E\ 3) (0)t\, 



where the following notations are used: 



(±1 



(P.*) 



-ifj 



(±i) 



• exp 



(-ib (±1 k) smid^H), 



(3.10) 



(3.11) 



(o) ( P ,t)= E 



In these formulas g( +1 ' — g, g( 1 ' = g* , 



;in(d (,T) t) +cos(d (CT) i) 



<T=±1 

x exp(-i6 (ff) t) 
-cxp[-rf (0) (p + A p )t] 



= [t(p±2A p )-i(p)]/2, 



(3.12) 



and 



6 (±i) = a (±i) + t(o)(p + Ap ), 



It is seen from Eqs. (|3.10| ) - (|3.12|) that an atom with an initial momentum component (along A p ) p ~ A p will 
change it to p ~ — A p at a time r„ (the time of the nir pulse) 



2 lsl 



(3.13) 



This transition is velocity-selective with the maximum efficiency determined by the Bragg resonance condition ( |2.S| ). 
The width of the peak in momentum distribution (the interval from the maximum to the first minimum) depends on 
the interaction time and for t = tq is 



6p = V3U\g\M/A p . 



(3.14) 



For a given A p it decreases with \g\. Therefore one should use a large detuning and small Rabi frequencies to get a 
narrower peak. 



IV. RAMAN EXCITATION 



Consider a three-level atom placed in the field of two plane electromagnetic waves ( p.l| ) with different frequencies 
(u>i — u)2 ~ whfs)- As before, the re gim e of large detunings is expected (|A{ 1-2 }| S> 7), which allows to neglect 



spontaneous emission and to use Eq. (3.2) for finding the density matrix evolution provided that G? oa '(pi, p^, t) is 
interpreted as the Green function of the three-component Shrodinger equation describing atomic dynamics during the 
|1) <-> 1 3), 1 2) «-> 1 3) transitions. 

If the detunings are also large enough in comparison with the maximal spatial shifts of transition frequencies and 
the laser intensities are far below saturation, i.e., the conditions 

|A 2 |»|fi 1 |,|fi 2 |,|f 3 -f 2 |L/?i, (4.1) 



G 



and are satisfied, the excited state of the atom can be eliminated adiabatically in analogy with Eq. (3.5). Making 
these and the rotating wave approximations one can reduce the Shrodinger equation so that it will involve only the 
wave functions ip±(p,t) and ip2(p,t) of atomic motion in the states |1) and |2) correspondingly. When rewritten in 
terms of closed family wave functions 0,|l8| (m = 1,2) 



^m{p,t) = exp[(i/7i).E|3>(0)£- iuj m t]ip m (p - hk m ,t), 



(4.2) 



this equation takes the form 



i-^ipi(p,t) = ti(p)ipi(p,t) +92ll>2 (p,t), 



(4.3a) 



d 

i-gifaiP^) = [t 2 (p) - ih ■ V]V>2(p,*) +3i^i(p,t), 



where g\ = fil^i/Ai, 92 — ^1^2/^2, and 

t m (p) =t(p- Kk m ) + A, 



I Or. 



(4.3b) 



(4.4) 



A. Quadrature solution 



Below the exact solution of the set of equations ( 4.3a ) - ( 4.3b ) is presented. First, the Laplace transformation is 
taken (m — 1,2) 



V>m(p,A) = / dter xt ^ m (p,t) 



(4.5) 



with the initial conditions ip m (p,t = 0) = ip m (p). Second, a Cartesian coordinate system in momentum space is 
introduced p = {p x ,Py,Pz) with the direction of z-axis chosen opposite to the vector f2. Then the equations for the 
Laplace transforms read as follows (/2 = | | ) : 



[ti(p) - iA]?/>i(p, A) + 32^2 (p, A) = -#i(p), 



(4.6a) 



d 

h{p) +2/2-5 »A 



fate, + gii>i(p, a) = -# 2 (p). 



Expressing ^(p, A) via ip 2 (p,ty from the Eq. (4.6a) 

-92ip2(p,X) ~ #i(p) 



V>i(p,A) 



h(p) - iA 



(4.6b) 



(4.7) 



one can simplify Eq. (4.6b) so that it will contain the only unknown function ^(p, A) and will allow an easy solution 
after imposing an appropriate boundary condition. This condition may be found if we take into account that an atom 
initially having a finite z-component of momentum is not able to reach infinitely large positive p z at any time because 
of the force fi acting in opposite direction, i.e., it is necessary to put 



^(p, A)| p ^oo = 0. 



(4.8) 



In such a way one gets 



1 f°° 
MP, A) = T I dp' z 
J2 J Pz 



Mp') 



9i4>i{p') 



h(p')-iX 



xe j^iHPiP'A) -iKvz -p' x )] 



(4.9) 



where p x and p y are fixed, i.e., p' = {p x ,Py,P z ), and 
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h(p,p',X) 



dp'! 



P'z 



Hp") 



9i 92 



h(p")-iX 



(4.10) 



The inverse Laplace transformation of expressions (4.7) (4.9) with the Mellin formula produces the desired quadra- 
ture solution (m = 1,2) 



ip m (p,t) = 2ni 



e+ioc 



dXe xt ip m (p,X), e>0. 



(4.11) 



B. Approximate formulas 



Although the integral in Eq. ( 4.11 ) cannot be calculated explicitly, it becomes possible to tabulate it after employing 
the following approximation. Note that by the retardation theorem the value of ^{1.2} (P? t) 1S nonzero only if t + 
(Pz — Pz)//2 > 0> which has a clear physical interpretation: due to the action of force f% a state with the momentum 
p z may arise only from the states which span the interval p z < p' z < p z + f2t. At times considered here (~ time of 
the 7r pulse) this interval appears to be narrow in comparison with the typical momentum in the system / 2 i <C A„. 
Therefore one can expand the function h(p, p', A) as a power series in p z — p' z and retain only the linear term 



h(p,p', A) ~ (p z -p' z ) 



*a(p) 



t\{p) - iX 

The validity of such an approximation is restricted by the condition 



t « JnhM/(f 2 A p ) 



(4.12) 



(4.13) 



preserving the phase of integrand in Eq. (4.9) against considerable variation due to omitted terms in h(p,p' , X). In 



the case of Na it requires the interaction time to be much less than 3.8 x 10 c. 

After changing the order of integration over dp' z and dX, subsequent calculations are obvious and lead to the 
following expressions for the Green function components (m,n = 1,2): 



G mn (p, p',t) = D mn (p + hk m ,p' + hk n ,t) 
xd±[p + fi(k m -k„) -p'] 
x exp[iuj m t - (i/h)E ]3) (0)t], 

which determine the evolution of wave functions ipi± t 2}(jp, t). Here 5j_(p) = S(p x )S(p z ) and 

Du(p,p',t) = exp[— iti(p)t]6(p z -p' z ) 

-V9i92V(P,p' ,t)X(+)(Pz -p' z ,t), 



(4.14) 



(4.15a) 



D 2 2(p,p , ,t)=exp[-it 2 (p)t]5(p t + f2t-j/ j ,) 

-v / 5i52'7(p,p',i)X(-)(Pz ~p' z ,t), 



(4.15b) 



Di 2 (p,p',t) = -ig2V(p,p' ,t)X(o)(p* -p' z ,t), 



(4.15c) 



D 2 i(p,p',t) = -igiv(p,p',t)x( ){Pz-p' z ,t). 



(4.15d) 



The values T)(p, p',t) and X(k)(Pz,t) (k — 0, ±) used in these equations contain the theta-function 0Q and the Bessel 
functions of the first kind J n Q (n = 0, 1) 



V(p,p',t) = exp| — (p z -_p^)[t 2 (p) - h(p)] - ih(p)t 
x^9(p z + f 2 t-p z )9(p' z - Pz ), 

J2 
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X(0){Pz,t) = J UV9l52ll - f 2 (p z ,t)]j 
1±/(P*,*) 



X(±)(Pz)*) 



X Jl {W9M1- P(Pz,t)}} , 



where /(p,,*) = 1 + 2p z /(f 2 t). 



C. Velocity selectivity 



To involve all atoms into the DS filling one must use a set of Raman pulses during every cooling cycle adjusting their 
number and the difference of beam frequencies for each pulse repetition in accordance with the velocity selectivity of 
discussed transition |1) — > \2), which appears to depend, besides other factors, on the size (~ 2L) of the atomic sample 
or, more precisely, on the width of the initial coordinate distribution. Indeed, any finite coordinate distribution, if 
fitted with the Gaussian profile oc exp(— r 2 /L 2 ), gives rise to a factor exp[— L 2 (p[ — p' 2 ) 2 / (4?i 2 )] in the density matrix 
Pn(p'ii P2: 0)- When (Lf 2 t/Ti) 2 C 1 (a thin atomic sample), this exponential factor does not influence the integration 
over dp'^ 2 j in Eq. (3.2), and for a flat initial atomic momentum distribution one gets standard formulas describing 



velocity-selective transitions in a three-level system without an external potential (see, e.g., [H]). In this case, after 
7r-pulse time 

to = 71-/ (2y/g^E) (4.16) 
the width of selection (half- width at 1/e of maximum) becomes 

5&P ~ h^^M/A p . (4.17) 



Conversely, at large L (when one can put pi ~ p 2 ~ p) the damping of the integrand in Eq. (3.2) is more rapid than 
phase variation of the Green functions product oc exp{i[i2(p) — ti(p)](p' lz — p 2z )/f2\, an d the velocity selectivity, 
which at small L has been provided by the resonance denominator, disappears. Instead, if h/L <C f 2 t, in the closed 
family basis it is determined by the exponential factor ~ exp{— h 2 [ti(p) — t 2 (p) ] 2 / (L 2 f 2 )} because of negligible theta- 



function contributions to the integral. Such picture is in agreement with Eq. (2.7). Thus, the width of selection may 



be evaluated as [cf. Eq. (2.8)] 



6 ( Vp~MLf 2 /(2A p ). (4.18) 

To illustrate this point in a ID configuration (when ki |J, k 2 and A p fj, f 2 ) let us regard a population p 2 2 = 
P22{p,PiTq) of the state \2) as a function of the momentum component p = p z and the dimensionless parameter 
i9 = Lf 2 To/h. Figure || shows the corresponding dependence calculated for sodium atoms with the initial density 
matrix pn(pi,P2,0) = exp[— L 2 (pi — p 2 ) 2 / (Ah 2 )] after 7r-pulse time. The peak characterizing velocity selectivity of 
the Raman excitation decreases with increasing i9, and begins to widen from the point 1? = n where 5^'p = S^p. So 
the parameter 1} d iscrim inates between the two domains with different behaviors of the selection width: when ?9 <C 7r 
it is given by Eq. ( 4.17j ) and docs not depend on the size of the atomic sample, whereas at 1? ^> ir it broadens with L 



in accordance with Eq. (4.18). In any case however, the width of selection cannot be less than S^p. 



V. OPTICAL PUMPING 



In this section, the optical-pumping beam is considered to be a plane electromagnetic wave (3J) at m — 2, whose 



interaction with atomic sample is described b y m eans of a master equation for the density matrix in external potential 
(see, e.g., [0). In the closed family basis |ll]Jll|, where the density matrix is denoted <7 a &(pi, P2, t), a,b = 1, 2, 3, and 
employing the rotating wave approximation, the master equation may be rewritten as 

d 



i-K7CT ab (p 1 ,p 2 ,t) = ^ [Hac{Pl)Sbd - H^ b (p 2 )5 ac ] 



cd 

xa cd (pi,p 2 ,t) +K ab {a}. (5.1) 



9 



The Hamiltonian H ac (p) is non-Hermitian because it governs the damping due to spontaneous decay of the excited 
state 



H a b(p) = [f(p - S a2 Kk 2 ) + A Q - i(5 a2 f 2 + 6 a3 i 3 ) ■ V] S a b 

-i("f/2)S a3 S b3 - fl 2 S a3 S b2 . 



(5.2) 



The last term in Eq. ( 5T ) is responsible for the return of the atom to the ground state via spontaneous emission 



T^abW] =5 a \5bil j dn$(n) 

xo-33(Pi - fik 3 ih,p 2 - hk 31 n,t), 



(5.3) 



where the function ^(n) determines the relative probability of emitting a photon with the wave number k 3 i — [E\ 3 \ (0)— 
E\i)(0)]/(hc) in the n direction. Below, the spherical symmetry approximation is adopted for <5>(n): $(n) ~ l/(4ir). 

Note that independent of the shape of excited-state distribution a 33 , the profile of atoms repumped into the state 1 1) 
is represented by a smooth functional 7?. a b[cr] whi ch undergoes sufficient variation only when its momentum argument 
changes by ~ hk 3 \ as it is evident from Eq. (5^3). A powerful simplification of the master equation may be achieved 
if we take into account that neither of the forces f 2 and £3 affect this functional significantly when the spatial shift of 
transition frequency is small enough 



max(|A 2 |, 7 /2)»|f 3 -f 2 |L/a. 



(5.4) 



Indeed, the last condition protects the excited-level population from the influence of the external potential, and the 
effect of forces reduces to a momentum kick dpk ~ f 2 T p + f 3 Td received by an atom, where r p and Td — I/7 are the 
lifetimes in the states |2) and |3) respectively. Evaluating t p from the optical pump rate (see, e.g., PJT^]) as 



7^+4A| 
4|f> 2 | 2 7 : 



(5.5) 



one can find that for a realistic set of parameters the value of the kick appears to be too small to induce a noticeable 
variation in the distribution of repumped atoms: |<5pfc| <C Tik 3 i. For example, in the case of sodium \Spk\ — 5xl0~ 6 hk 3 i 
if IO2I = O.I7, and A 2 = 0. Therefore the potential energy term containing forces f 2 and f3 can be omitted in the 
Hamiltonian (5.2). 

After such a simplification the master equation admits an analytical solution, which one can obtain, e.g., with the 
help of Laplace transform. As a result, the only component of the original-basis density matrix, nonvanishing at 



t t p , takes the form 



Pn(Pi,P2,i) = /Oii(Pi,P2,0)exp(-iA t) - 



47T 



dh 



E E 

a=±l/3=±l 



(-1) 



a+0 



exp(-iA i) - exp{-it[A( a) (pi) - A^p^)]} 



In this equation p'^ 



2} 



A( Q )(pi) 
P{i,2} - ^ 31 n, and 



^(/3)( p 2) _ ^0 



t(Pl) ~ i (P2), 



-P22(p'i + ?ik 2 ,p 2 + ftk 2 ,0). 



(5.6) 



D(p) = [S^(p)f+A\n 2 



1 r 



\(±1) 



where 



S (±) (p) = t(p - Kk 2 ) + A 2 ± [t(p) - » 7 /2] 



The above solution is valid when the spontaneous decay rates into states other than |1) are negligible. Due to the 
specific choice of |3) as the working excited state all such decay rates turn out to be Tyfjp <C 1 times less than 7. 



10 



Nevertheless, in real atomic system some fraction of atoms will accumulate in unwanted states. To return them back 
into the ground level one should include additional laser beams in the considered scheme, as it is done, e.g., in the 
coherent optical pumping pp[ |. 

Another complication may arise from non-resonance excitations out of the ground level which we do not include in 
the treatment for both the second and third steps of cooling. In effect, the laser light driving the |2) — * |3) transition 
introduces a detuning ~ ^hfs and a Rabi frequency ~ with respect to the transition |1) — * |3). However, if 

we restrict f^2 by the condition 

^2<^HFS??HF7 (5-7) 

the excited-state population [oc ^/(whfs^hf) 2 ] appears to be small. 



VI. NUMERICAL SIMULATION 



In the following we present one-dimensional results obtained for Na assuming that all vectors have only z- 
components, i.e., lie on the same axis with the gravitational force, and the laser beams with ki and k2 are counter- 
propagating. An initial distribution of ground-state atoms is considered to be Gaussian 



Pii(pi,P2,0) 



L 



ft7rV2 J\ + 2L?ol 



: exp 



Ah 2 



(pi +P2) 2 
1 + 2L 2 ct2 



+ (pi ~ P2) 5 



(6.1) 



where p — p z , and o~ p stands for momentum dispersion (in units of K). Since in our scheme we imply that an atomic 
sample precooled to the recoil limit is used, it is reasonable to take the wave number of laser light k = 1.07 x 10 5 cm -1 
as a n input for <r p . Note that for the considered laser-beams geometry A p /h~ = k. In order to satisfy the condition 
(5.7) we also take the parameter 77HF = 2.5 x 10~ 2 , which corresponds to Bq = 10 3 G. 

In the numerical simulation of a cooling cycle each first-step pulse was followed by five repetitions of a set involving 
seven Raman and one optical pumping pulses. A map ping b etwe en the density ma trices at the beginning and the 
end of this sequence of pulses is given by the formulas (Q, (|3.1C|), (|4.l4) , a nd jO ). Both first and second steps of 

is. (fug) and (|T§c( 



cooling continued during r -times defined according to the Eqs. ( |3.13 ) and ( 4.16| j correspondingly. The duration of 
the optical pumping pulse was taken to be 10r p to provide a complete depopulation of the |2) state. The remaining 
parameters were chosen as follows. For the first step: Oi = Q2 = O.O87, and Ai = — 167. For the second step: 
ill = OT7, J7 2 = 0.047, and all seven Raman pulses were detuned to the red so that the sum + A 2 = —327 
remained constant while the difference Ai — A2 was increased by -135, 118, 372, 625, 880, 1135, and 1393 kHz. Such 
a choice of detunings was tailored both to span the momentum interval < p < 3A p and to minimize the losses of 
DS population due to parasitic excitation by sidelobes in the frequency spectrum of Raman transitions at ?9 < it (see 
Fig. H). For the third step we put Q2 = O.I7, and A2 = 0. The initial size of atomic sample was taken L = 1 cm. 
However, for the given set of Raman light parameters this, or indeed any sma ller, value of L leads to the inequality 
■d < tt. It means that the width of velocity selection is determined by Eq. ( 4.17 ) and does not depend on L. Therefore 
our results remain correct for all L < 1 cm. 

Figure [| shows the initial momentum distribution p\\ (p, p, 0) and the formation of a DS peak during two first cooling 
cycles including intermediate stages when the position of this peak is alternated. Although each |1) — > |1) transition 
captures atoms in a rather wide momentum interval ~ 2Sp w 0.28A p , the width of the DS peak (at half-maximum) 
decreases rapidly with the number of applied cycles because of a pronounced maximum in the transition rate profile. 
After 10 cycles the decrease slows down and approaches at 0.005A p by the end of cooling, as may be seen from Fig. 
^(a). At the same time, the peak height growth is far from saturation, as depicted in Fig. ^(b). After 100 cooling 
cycles the height is more than 98 times the initial distribution maximum. 

The fraction of cold atoms in the interval — A p — Sp < p < — A p + Sp depends on the difference between the feeding 
rate due to optical pumping and losses during DS transfer, which arise along with the reduction of the peak width. 
Figure |J(c) shows that about 65 % of all atoms collect there by the end of cooling. 

When separated from the background by the final first-step transition which transfers the aforementioned interval 
to positive momentum half-axis, the DS peak acquires a shape represented in the Fig. |5|. As a result, the effective 
temperature calculated as a mean kinetic energy of the atoms distributed within the domain A p — dp < p < A p + Sp 
reaches 0.4 nK or 0.015T G . 

A very special sequence of laser pulses considered in our numerical simulation was designed both to minimize the 
volume of computations and to demonstrate the possibility of cooling below the gravitational limit with a noticeable 
efficiency. However, this sequence is not the best from the point of view of practical application because it includes 
Raman pulses in the regime $ < n, in which excitation-spectrum sidelobes can destroy DS unless one correctly adjusts 
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detunings. For experimental purposes the alternative choice $ > tt may be much more attractive insofar as it does 
not lead to parasitic excitations, as illustrated in Fig. Unfortunately, in this case each cooling cycle must contain a 
considerable number of steps because of the moderate excitation rate of Raman transition, and a numerical simulation 
does not seem feasible. 



VII. CONCLUSIONS 



In this paper we have studied a modification of Raman cooling method, in which the atomic internal ground 
state possesses a translational invariance due to compensation of gravity by the Stern-Gerlach effect. Our scheme 
is based on creating a dark state which is cyclically moved in momentum space with additional velocity-selective 
two-photon transitions so that it is kept unreachable for Raman excitation-repumping pulse sequences at all times. 
The consideration has been restricted to dilute atomic samples, i.e. we have not included any many-atom interactions 
plf . I n the approximation in which the DS losses induced by atomic collisions, non-resonance excitations, and photon 
scattering are neglected, our one-dimensional numerical computations have shown that a considerable fraction of 
all particles can be cooled to the temperature below the gravitational limit. Furthermore, this temperature can be 
decreased significantly for laser pulses which provide a smaller width of the velocity selection during DS transfer. 
Though we have given a numerical example only for a one-dimensional model, our theoretical investigation of the 
suggested scheme is also valid for the two- and three-dimensional cases. 

In a bosonic system, where the losses in DS population can be compensated by the quantum-statistical enhancement 
of feeding rate, our scheme will be appropriate for dense samples as well. Thus, it may be used in creation of a coherent 
atomic-beam generator |l^,^2|. An easy tunable wavelength will be one of the advantages of such a device, because 
the momentum of a cooled atom is readily defined by the geometry of laser beams. 



[1] P. Lett, R. Watts, C. Westbrook, W. Phillips, P. Gould, and H. Metcalf, Phys. Rev. Lett. 61, 169 (1988). 

[2] A. Aspect, E. Arimondo, R. Kaiser, N. Vansteenkiste, and C. Cohen- Tannoudji, Phys. Rev. Lett. 61, 826 (1988); J. Lawall, 

S. Kulin, B. Saubamea, N. Bigelow, M. Leduc, and C. Cohen- Tannoudji, ibid. 75, 4194 (1995). 
[3] M. Kasevich and S. Chu, Phys. Rev. Lett. 69, 1741 (1992); N. Davidson, H. J. Lee, M. Kasevich, and S. Chu, ibid. 72, 

3158 (1994); H. J. Lee, C. S. Adams, M. Kasevich, and S. Chu, ibid. 76, 2658 (1996). 
[4] R. Dum, Phys. Rev. A 54, 3299 (1996). 
[5] R. Dum and M. Ol'shanii, Phys. Rev. A 55, 1217 (1997). 
[6] T. Pellizzari, P. Marte, and P. Zoller, Phys. Rev. A 52, 4709 (1995). 
[7] G. Morigi, J. I. Cirac, K. Ellinger, and P. Zoller, Phys. Rev. A 57, 2209 (1998). 

[8] J. Reichel, F. Bardou, M. Ben Dahan, E. Peik, S. Rand, C. Salomon, and C. Cohen- Tannoudji, Phys. Rev. Lett. 75, 4575 

(1995) . 

[9] P. E. Moskowitz, P. L. Gould, S. R. Atlas, and D. E. Pritchard, Phys. Rev. Lett. 51, 370 (1983). 
[10] P. J. Martin, B. C. Oldaker, A. N. Miklich, and D. E. Pritchard, Phys. Rev. Lett. 60, 515 (1988). 
[11] A. V. Soroko, J. Phys. B. 30, 5621 (1997). 

[12] A. P. Kazantsev, G. A. Ryabenko, G. I. Surdutovich, and V. P. Yakovlev, Phys. Rep. 129, 75 (1985). 
[13] W. Zhang and D. F. Walls, Phys. Rev. A 52, 4696 (1995). 

[14] E. A. Korsunsky, D. V. Kosachiov, B. G. Matisov, and Yu. V. Rozhdestvensky, JETP 76, 210 (1993). 

[15] W. Zhang and D. F. Walls, Phys. Rev. A 49, 3799 (1994). 

[16] A. M. Guzman, M. Moore, and P. Meystre, Phys. Rev. A 53, 977 (1996). 

[17] V. Minogin and V. Letokhov, Laser Light Pressure on Atoms (Gordon and Breach, New York, 1986). 

[18] A. Aspect, E. Arimondo, R. Kaiser, N. Vansteenkiste, and C. Cohen-Tannoudji, J. Opt. Soc. Am. B 6, 2112 (1989). 

[19] B. Matisov, V. Gordienko, E. Korsunsky, and L. Windholz, JETP 80, 386 (1995). 

[20] E. A. Korsunsky, W. Maichen, and L. Windholz, Phys. Rev. A 56, 3908 (1997). 

[21] M. Lewenstein, L. You, J. Cooper, and K. Burnett, Phys. Rev. A 50, 2207 (1994). 

[22] R. J. C. Spreeuw, T. Pfau, U. Janicke, and M. Wilkens, Europhys. Lett. 32, 469 (1995); H. M. Wiseman and M. J. Collett, 
Phys. Lett. A 202, 246 (1995); M. Holland, K. Burnett, C. Gardiner, J. I. Cirac, and P. Zoller, Phys. Rev. A 54, R1757 

(1996) ; G. M. Moy, J. J. Hope, and C. M. Savage, ibid. 55, 3631 (1997). 



12 




FIG. 1. Energy level diagram and profile of ground state atomic distribution pn as a function of p, the momentum 
component along the direction of vector A p : (a) before DS transfer with the two-photon |1) — > |1) transition; (b) before the 
Raman excitation |1) — + |2) cycle; and (c) before the optical repumping pulse tuned to the |2) — > |3) transition. The curly 
braces denote momentum intervals involved in each of these processes. 
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FIG. 2. Velocity selectivity of Raman excitation in ID case for Na represented as the dependence of P22, the state \2) 
population, on the momentum component p — p z and parameter 1? = Lf2To/h. Detunings Ai = A2 = — 167 and Rabi 
frequencies fl\ — O.I67, SI2 = O.OI67. 
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FIG. 3. One-dimensional ground-state distribution p\\ for Na as a function of atomic momentum p and the number n of 
cooling cycles normalized to 1 on the scale p/A p . Half-integer values of n correspond to the beginning of the second stage of 
each cooling cycle which starts when laser beams reverse. The curve with n = gives the initial distribution. The highest 
peaks of the function represent the DS. 
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FIG. 4. Parameters of the DS peak vs the number of cooling cycles, (a) Peak width at the half-maximum, (b) peak height, 
and (c) fraction of cold atoms around DS in the interval ±<5p ~ ±0.14A P , which spreads between the first minima in momentum 
dependence of the |1) — > |1) transition rate. 
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FIG. 5. Final ground-state distribution pu of atoms with positive momenta as a function of p. The full-scale distribution 
is normalized as in the Fig. ^| whereas the plotted part contains ~ 59% of all particles. 
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